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Abstract 


Single-molecule discrimination among amino acids is crucial to the realization of next-generation protein 
sequencing. Owing to the heterogeneous charge and subtle volume difference of underivatized amino acids, it 
remains a challenge for single-molecule techniques to recognize each of them. Here, we report the direct detection 
of twenty proteinogenic amino acids using a copper(II)-functionalized MspA nanopore. The binding sites for 
copper(II) ion are constructed by introducing histidine mutation (N91H) to M2MspA protein. With copper ion 
binding to histidine residues, amino acids can reversibly coordinate the copper-histidine complex, generating 
well-defined current signals. Using this strategy, all twenty amino acids can be detected. Assisted by a machine 
learning algorithm, we can identify 100% of signals with 70.2% accuracy or 60% of signals with 93.4% accuracy in 
the validation set. In successively addition experiment, each amino acid in a mixture of 10 amino acids can be 
identified precisely. Furthermore, we use carboxypeptidase Al to partly release the C-terminal amino acids of 
peptides with different lengths (9, 10 and 22 residues). The hydrolysates of peptides can be identified and 
distinguished. These results demonstrate the feasibility of this system for amino acids detection and peptide 


identification, shedding new lights on the development of single-molecule protein sequencing. 


Introduction 


Amino acids are building blocks of proteins, raw materials for biosynthesis, playing fundamental roles in various 


physiological and pathophysiological processes such as epigenetic regulation and tumor metabolism!~. Therefore, 
it is crucial to detect and identify amino acids with a higher spatiotemporal resolution, which has recently aroused 
great interest for researchers, especially in the field of single-molecule protein sequencing>*. Due to alternative 
RNA splicing and post-translational modification, the resulting proteoforms are highly complicated and contain 
deeper-level information that cannot be accessed directly from transcriptome’. In addition, there is no existing 
method to amplify proteins similar to DNA amplification. Consequently, it is difficult for mass spectrometry-based 


10,11 


methods to identify low-abundance proteins from proteome’’’. To address these problems, single-molecule 


methods that can distinguish the twenty proteinogenic amino acids should be developed for protein sequencing. 


For fluorophore-based methods, specific amino acids like cysteine and lysine can be selectively modified by 
fluorescent molecules. Then, by sequentially degrading peptide using Edman chemistry, or direct imaging using 
single-molecule FRET, the relative position of labeled amino acids can be deduced from the fluorescent signals!?""*. 
Additionally, fluorophore-labeled N-terminal amino acid recognizers have been engineered to bind specific amino 


acids reversibly!5!6 


. The repetitive signals of the same amino acid can greatly improve the accuracy of 
identification!’. Although these methods have high throughput and good reliability, it is difficult for chemists to 
label twenty amino acids specifically. For label-free methods, techniques such as tunneling current 


measurement!*!? and molecular junctions’ enable rapid and precise detection of up to 12 amino acids. 


Given the nanopore technique has demonstrated its superiority in single-molecule DNA sequencing, it has been 
also considered an ideal candidate for amino acid detection and protein sequencing>!””. Studies have shown that 
peptides with different properties can be directly detected and distinguished, such as molecular weight”>4, 


27.28 and single-amino acid mutation”’. To further analyze the amino acid 


length*>”°, post-translational modification 
sequence of the peptide, the translocation of the peptide must be well controlled to generate sequence-dependent 
signals. Protein unfoldase ClpX was used to unfold and drive protein through a nanopore, and different segments of 
protein could be discerned*®. Moreover, the ratcheting motion of DNA-peptide conjugation through nanopore was 


achieved using DNA helicase or polymerase, generating clear sequence-dependent signals*!*?. 


Single-file 
translocation of linearized protein can be facilitated by engineering electroosmotic force***>. Unfortunately, it is 
challenging to deconvolute the signal contributed by 5-6 amino acids, as twenty types of amino acids generate more 
combinations than four types of nucleotides. Therefore, analysis of single amino acid translocation events could 
provide valuable information. Underivatized amino acids can be detected using copper ion-modified a-hemolysin, 


MoS>2 nanopore**?’ 


. Furthermore, taking advantage of the pore structure, the aerolysin nanopore is capable of 
distinguishing 13 out of 20 amino acids only with polyarginine carrier*®. Biosensors that can discriminate all twenty 
proteinogenic amino acids directly with high sensitivity and specificity at the single-molecule level are required for 


biosensing and protein sequencing. 


Here, we report the direct identification of twenty proteinogenic amino acids using a copper(II)-functionalized 
MspA nanopore. Benefiting from conical pore geometry, MspA nanopore has proved to be an ideal choice for 
sensing ions and small molecules*®*°. We constructed binding sites for copper ions by introducing histidine 
mutations in the restriction region of the pore lumen. With copper ion binding to histidine residues, the reversible 
coordination between amino acid and copper-histidine complex could generate well-defined current signals. Next, a 
machine learning-based classifier was trained for the identification of twenty amino acids. Furthermore, three types 
of peptides (a-Bag Cell Peptide (1-9), adrenocor ticotropic hormore (ACTH, 18-39) and Angiotensin I) were partly 
hydrolyzed from C-terminal using carboxypeptidase Al. The composition of their hydrolysates could be identified 
accurately. These results suggest the great potential of this system for amino acid detection and peptide 


identification, paving the way for next generation protein sequencing. 
Results 


Sensing of twenty amino acids using copper(II)-MspA 

In a typical experiment, amino acids and copper(II) chloride are added to the cis (grounded) and trans chamber 
respectively (Fig. la). The binding sites for copper ions are located at the constriction region of MspA nanopore 
(Fig. 1b). For each of the eight subunits, the 91* asparagine is substituted by histidine to create a copper-binding 
structure which is similar to histidine brace motif*!. We suppose that one 90" asparagine residue and two adjacent 
91%* histidine residues could reversibly coordinate a single copper ion and then amino acid (Fig. 1c). The 
corresponding three binding states could be observed from stepwise current change (Fig. 1d). According to this 
supposition, there are at least four binding sites for copper ions, and the reversible binding of multiple copper ions 
were also observed}. However, such stochastic binding events interfered with the precise assay of subsequent 
amino acid binding. To keep the current baseline at one constant level (Io), excess copper ions at a final 
concentration of 200 uM were added to saturate the binding sites during most of the measuring time (c.a. 87.8 + 
3.1%) (Supplementary Table 1). 


All twenty amino acids can be detected and produced clear signals with high reproducibility (Fig. le). Commonly, 
blockade (|Ii-Io|/Io) and dwell time (At) are analyzed to characterize the signals. The signal blockades for each type 
of amino acid exhibit a unimodal distribution (Fig. 2a). Most of them can be well distinguished from each other. 
The mean blockade and its standard deviation of each type of amino acid were calculated from the mean value of 
the Gaussian fit. The blockades are shown to have a good positive correlation with the volume of amino acids (Fig. 
2b). When amino acids with charged side group and proline (P) are excluded, the coefficient of determination of 
linear fitting reaches up to 0.92 (Supplementary Fig. 1), indicating that the generation process of current blockade 
for these amino acids obeys the classical volume exclusion model. While for amino acids with charged side groups, 
the volume exclusion model is no longer applicable**. The blockade of aspartic acid, glutamate acid and histidine 
(D, E and H) is bigger than expected. It could be attributed to the possible interaction between their side chain and 
copper-histidine complex. For lysine and arginine (K and R), the electrical repulsion between their amine groups 
and copper ions was expected to lower the current blockade. Due to the strong interaction between copper ions and 
the sulfhydryl group of cysteine (C), the binding of copper ions to histidine residues could be extremely unstable 
(Supplementary Fig. 2). The fluctuation of open pore current made it hard to determine the baseline current Io and 
also shortened its duration. Therefore, few signals of cysteine were extracted, causing a high standard deviation of 


the mean blockade. 


In terms of signal frequency, there are remarkable differences among different amino acids (Fig. 2c, e). Among 
them, the signal frequency of P is the lowest because of its unique structure, which can hardly interact with copper 
ion. The mean signal frequency of polar amino acids is significantly higher than that of nonpolar amino acids. For 
amino acids with charged side chains, signal frequencies of K and R are significantly lower than that of negatively 
charged D and E. One of the reasons may be that K and R can be driven away from the nanopore by electrophoresis 
force. The isoelectric point of H (7.59) is quite close to 7.5. And the interaction between H and copper ion is very 
strong, so its signal frequency is even higher than that of D and E. Therefore, the charge distribution of amino acids 
and their interaction with copper ion both contribute to the signal frequency. To demonstrate that the a-amine group 
and a-carboxyl group of amino acids are essential for coordination, we synthesized two dipeptides, EF and yEF 
(peptide bond is formed using the y-carboxyl group of glutamate). The translocation of dipeptide yEF generated 


distinguishable signals with larger blockade (37.8 + 0.6% vs 15.3 + 0.2%) and longer dwell time (2.521 + 0.149 ms 
vs 0.554 + 0.014 ms) than dipeptide EF (Supplementary Fig. 3). Since all these amino acids bind to a copper ion 
with their a-amine group and a-carboxyl group, the mean dwell times of signals are in the same order of magnitude, 
ranging from 1.87 ms (G) to 6.86 ms (W) (Fig. 2d). Additionally, the acetylated leucine and amidated leucine did 
not generate characteristic signals, which provided direct evidence for our chemical model (Supplementary Fig. 4). 


Amino acid identification by machine learning 

Multiple amino acids can bind to nanopore at the same time, resulting in superimposed multi-level signals 
(Supplementary Fig. 2). These signals cannot be identified simply by the two parameters of the blockade and dwell 
time. To improve the accuracy of amino acid identification, we trained a machine learning-based classifier (Fig. 3a). 
We firstly normalized the classified amino acid signals through dividing current amplitude of signals by their 
baseline current (Io). The distribution of current density of each normalized signal was divided into 1000 equally 
sized intervals from 0 to 1. Features X0001-X1000 were then calculated from data of their corresponding interval, 
representing the density of data points within each interval. Then, the mean blockade, dwell time, standard 
deviation and features X0001-X1000 of normalized signals were used as input features to train the classifier using a 
machine learning algorithm. Next, we assessed the performance of six different classification algorithms. Results 
showed that the random forest (RF) algorithm worked best, and the corresponding receiving operator characteristic 
(ROC) analysis revealed an area under the curve (AUC) of 0.9708. Features ranging from X150 to X178 have 
larger importance values compared with others (Fig. 3b). These features were generated from signal points within 
the range of level 1 blockade of amino acids, indicating that the part of the signal generated from the single amino 
acid binding event was still the most important part for signal identification. 


To minimize the influence of background signals and further improve the accuracy of amino acid identification, we 
filtered signals by applying cutoff value to dwell time. Signals with dwell time lower than 1 ms, 2 ms, 3 ms, 4 ms, 5 
ms, and 6 ms were filtered in RF1-6 models, respectively. We found that with the increasing cutoff value, the 
performance of models improved, despite the reduction in the number of available signals (Supplementary Fig. 6a, 
b, c). The results indicate that signals with higher dwell time are more likely to be correctly identified. The AUCs 
of ROC in the test set and validation set increased from 0.9837 and 0.9708 to 0.9936 and 0.9838 respectively 
(Supplementary Fig. 6d, e). In addition, the sensitivity, specificity, precision, recall and F1 score of random forest 
classification models improved significantly (Supplementary Fig. 7a), and the prediction probability of corrected 
signals is also enriched around 1 (Supplementary Fig. 7b). 


The classification accuracy of the test set ranged from 75.64% to 87.66% with 100% signal recovery. When a 
prediction probability cutoff (>0.7) was used to filter out unclear results, the average classification accuracy 
increases up to 95.64% with 63.85% signal recovery (Supplementary Fig. 6f). In the validation set, the accuracy 
ranged from 70.2% to 83.2% with 100% signal recovery. When only signals with a prediction probability greater 
than 0.7 were considered, the average classification accuracy increases up to 93.4% with 60% signal recovery (Fig. 
3c, Supplementary Fig. 7b). These results indicate that there is no over-fitting in these models. In model RF4, the 
accuracy reached up to 95% with all amino acids showing good differentiation except C, N, and T (Fig. 3d). 
Furthermore, in successively addition experiment, the results showed that each amino acid in a mixture of 10 
proteinogenic amino acids and S-carboxymethyl-L-cysteine can be identified precisely (Supplementary Fig. 5). 


Discrimination of peptides by recognizing their hydrolysates 
As it is challenging to sequence a linear peptide directly, the detection of individual amino acids cleaved from a 


peptide may offer an alternative. We next tested the feasibility of this system for the identification of peptide 
hydrolysate. Carboxypeptidase Al was used to sequentially release single amino acid from C-terminus of peptides. 
The hydrolysis reaction stops when either of three amino acids (K, R and P) becomes the first amino acid at 
C-terminus (Fig. 4a). In this way, it can not only avoid the detection of K, R and P, but also produces limited types 
of amino acids, reducing the complexity for amino acid identification. The results showed that the hydrolysis of 
peptides could be monitored in real-time, after mixing carboxypeptidase Al with peptide EAFNL directly in cis 
chamber (Fig. 4b). To make the hydrolysis reaction more sufficient, three types of peptides (a-Bag Cell Peptide 
(1-9), ACTH (18-39) and Angiotensin I) were hydrolyzed respectively with a higher concentration of 
carboxypeptidase Al at a temperature of 37 °C, and then added to cis chamber. Their hydrolysates were detected 
and identified, which were mostly consistent with the theoretical amino acid products (Fig. 4c). For all three types 
of peptides, the percentage of leucine signals is lower than other theoretical products (Fig. 4d), which was similar to 
the trend of signal frequency when amino acids were detected separately (Fig. 2c). For ACTH (18-39) and 
Angiotensin I, some of their hydrolysates were identified incorrectly, as F and L were identified as other amino 
acids with a close blockade. Although the hydrolysates of these peptides can be distinguished, the accuracy for 
identifying amino acids mixture still needs to be improved. 


Discussion 


In summary, this study enables the direct detection of twenty proteinogenic amino acids using a copper ion 
functionalized MspA nanopore. The interaction between o-amine group, a-carboxyl group of amino acid and 
copper-nanopore complex is considered the key to generating current blockade. Consequently, unnatural amino 
acids are supposed to be detected. A shortcoming of this method can be multiple binding sites for copper ions and 
amino acids due to the sequence homology of eight MspA subunits, which resulted in superimposed event signals 
(Supplementary Fig. 2). Therefore, MspA with two adjacent N91H-mutant subunits is required for the binding of 
single copper ion and amino acid, which helps to reduce signal complexity and improve sensing accuracy***>. In 
addition, we noticed the maximum blockade of amino acids only accounted for ~1/4 of open pore current. This 
could attribute to the large pore diameter and short translocation duration of amino acids. So a narrower pore with 
stronger interaction with amino acids should be rationally designed*®*’. It is observed that besides hydrophobic 
volume, the charge of amino acids has considerable influence on the current blockade. This needs to be further 
investigated combined with theoretical calculation and may provide a valuable reference for polypeptide 


sequencingt 4s 


. Compared with analysis using blockade and dwell time as parameters, the machine learning 
algorithm developed here allows the analysis of all data points of one signal, improving the accuracy for amino acid 
identification. We also demonstrated the C-terminal amino acids of the peptide could be partly released by 
carboxypeptidase Al and identified using this system. As single-molecule protein sequencing would revolutionize 
proteomics research and has valuable practical applications such as identifying neoantigens, this method needs to 
be further developed towards single-molecule resolution. It may be beneficial to construct a peptidase-nanopore 
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conjugation or even a nanopore with peptidase activity“, which could enable the hydrolysis and sequencing of a 


single peptide. 
Methods 


Protein preparation. M2MspA-N91H mutant was expressed and purified as described previously**!. Briefly, A 
gene of M2MspA with 91“ histidine mutation for each of eight subunits was cloned into pET28b vector. Then, the 
plasmid gene was heat-shock transformed into E. coli BL21 (DE3) competent cells. The cells were cultured in LB 


medium containing kanamycin (50 wug/mL) to an ODoeo of 0.8, and then 1 mM isopropyl 
B-D-1-thiogalactopyranoside (IPTG) was added. Afterward, cells were incubated at 15 °C for 12 h with 220 rpm 
shaking. The cells were harvested by centrifugation at 4000 rpm, 4 °C for 15 min and then re-resuspended. Cell 
disruption was performed by Sonication using an Ultrasonic cell disruption device. Then the supernatant was 
retained, and the target protein was furtherly purified using anion exchange column (Q-Sepharose) and size 
exclusion column (Superdex 200 16/90). 


Amino acids detection and peptide hydrolysis. Electrophysiology experiments were performed using a classical 
vertical lipid bilayer setup (Warner Instruments). A pair of Ag/AgCl electrodes were placed in the trans and cis 
(grounded) side of the chamber, which was filled with 1 mL of electrolyte solution (1 M KCl, 10 mM MOPS, pH 
7.5). Then, the planar lipid bilayer membrane was formed on the 150 um-diameter aperture by painting a thin film 
of 1,2-diphytanoyl-sn-glycero-3-phosphocholine (DPhPC) (Avanti Polar Lipids). A voltage of +300 mV was 
applied to induce nanopore insertion after adding the MspA protein (final concentration of 60-90 ng/mL) into the 
cis chamber. After a single nanopore insertion, CuCl2 solution was added into the trans chamber to a final 
concentration of 200 uM (20 uM in peptide hydrolysis experiments). High-purity L-amino acids (>98%) were 
dissolved in Milli-Q water away from light immediately before use. To collect more signals, amino acids were 
added to the cis chamber to a high final concentration of 100 uM (except 5 uM, 200 uM, and 2 uM for H, P and C, 
respectively). For peptide hydrolysis, the peptide was dissolved in Milli-Q water to a final concentration of 2 mM. 
8 uL peptide solution was mixed with 2 uL 3.3 U carboxypeptidase Al and reacted at 37 °C for 15 mins, then the 
mixture was added to the cis chamber. For real-time monitoring of peptide hydrolysis, peptide N’-EAFNL-C’ was 
added to cis chamber to a final concentration of 20 uM, followed by the addition of 10 uL 16.7 U carboxypeptidase 
Al. 


Electrophysiology recording. Single-channel current recordings were amplified using an Axopatch 200B amplifier 
(Molecular Devices) and filtered with a built-in four-pore low-pass Bessel filter at 2 kHz. Data were digitized by a 
Digidata 1550B converter (Molecular Devices) at a sampling rate of 100 kHz. All electrophysiology experiments 
were performed at room temperature (23 + 2 °C). 


Signal extraction for amino acid translocation event. Firstly, to reduce the noise of raw current recording, we 
calculated the optimal changepoints of current according to the mean and variance and polished the current 
recording using the mean current of each segmentd time range according to the identified changepoints. Then we 
extracted the translocation events from the polished signal based on the minimum blockade threshold value (0.1) 
against the baseline current. For all the extracted events, we calculated the blockade, dwell time, and standard 
deviation of signal current (SD) et al.. In addition, in order to better describe the characteristics of each signal, we 
uniformly extracted the density values of 1000 points from the density curve of the normalized current of each 
signal as the characteristic values of the signal. 


Raw signal filtering based on the similarity with background noise. For the raw signals of each independent 
experiment, we randomly selected 

The extracted raw signals was 

For the original signal of each independent experiment, we randomly selected the same number of noise signals 
from the corresponding blank control to calculate the Euclidean distance matrix. The eigenvalue of Euclidean 
distance is the predicted value of the machine learning model. 


Classification model training. We developed a machine learning (ML) algorithm to automatically predict the 

corresponding amino acid from the signal of a translocation event. The strategy was to utilize an algorithm to “learn” 
from the classified training data and build an optimum classification model to recognize unknown events. To train 

the model, the blockade, dwell time, SD value and the estimated density values of the normalized signal were 

calculated using R program to form a feature matrix (Fig. 3). For each amino acid, we randomly selected one of the 

experimental data as the independent validation set, and then randomly selected 80% of all the remaining signals as 

the training data set, and 20% as the test set. Model training was performed using the R package caret. A set of 

classifiers including random forest (RF), naive Bayes (NB), K nearest neighbours (KNNs), Bagged CART, 

AdaBoost.M1 (AdaBoost), and neural networks (NNet) classifiers were tested. To prevent the over-fitting of model 

training, 10-fold cross-validation was performed for each model to report the cross-validation accuracies. 

Considering that the dwell time of the signal reflects signal quality, and signal with dwell time lower than 1 ms may 

be generated from spontaneous gating of the nanopore, signals were then filtered in each model to improve the 

accuracy. Signal with dwell time lower than | ms, 2 ms, 3 ms, 4 ms, 5 ms, and 6 ms was filtered in model RF1, 2, 3, 
4, 5, and 6, respectively. Finally, the trained model was then used to predict unclassified events. 
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Figure 1. Experimental setup and principle of amino acid detection. a, Schematic illustration of experimental 


setup. Amino acids and copper ions were added to cis and trans chamber respectively. A voltage of +50 mV was 


applied during measurement. The N91H mutant sites of eight subunits are highlighted in orange. b, Bottom-view 
structure of M2MspA-N91H nanopore (predicted using SWISS-MODEL). The dotted box shows a binding site for 
copper ion. c, Proposed sensing mechanism. Two adjacent histidine residues and one 90" asparagine residue firstly 


coordinate a copper ion. The o-amine and a-carboxyl group of amino acid then coordinate the histidine-copper 


complex. d, A representative current trace shows the corresponding current change for three binding states in Fig. 


lc. e, A typical current signal of amino acid translocation event. 
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Figure 2. Characteristics of signals of twenty proteinogenic amino acids. a, The distribution of relative 
abundance of the blockade of amino acid signals. (n = 4278 (E), 4211 (D), 650 (K), 193(His1), 306(His2), 7166 (F), 
3934 (W), 2768 (Y), 3025 (I), 8004 (M), 3059 (R), 8131 (T), 8101 (S), 3750 (L), 857 (A), 1149 (G), 361 (P), 7873 
(Q), 9634 (N), 2119 (V), 616 (C)) b, Scatter plot of volume versus blockade of amino acids. For each amino acid, 
the blockade and its standard deviation was calculated using mean values of Gaussian fit from at least three 
independent experiments. Amino acids with charged side chain, nonpolar and polar amino acids are colored in 
green, purple and orange, respectively. Pearson correlation coefficient = 0.58, p = 0.0062. c, Boxplot of signal 
frequency of amino acids. Each dot represents data from an independent experiment. d, Mean dwell time of amino 
acid signals. e, Boxplot of signal frequency of four categories of amino acids. The signal frequency of polar amino 
acids is significantly higher than nonpolar amino acids. 
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Figure 3. Amino acid identification assisted by machine learning algorithm. a, Illustration of training process. 
First, classified level 1 (one amino acid binding) and level 2 (two same amino acids binding) signals of each type of 
amino acid were imported and normalized. Then the level 1 blockade, dwell time, standard deviation were 
extracted. Additionally, 1000 data points were extracted from the current density of each signal (from 0 to 1 with an 
interval of 0.001), named Feature X0001-X1000. Performance of models were tested including random forest (RF), 
naive Bayes (NB), neural networks (NNet), K nearest neighbour (KNN), Bagged CART and AdaBoost.M1 
(AdaBoost). Among these, random forest model was the best one with AUC of 0.990. A 10-fold cross-validation 
was used to prevent over-fitting. b, Feature importance generated from training of random forest model (RF) for L1 
signals of all 20 amino acides. The upper y axis represents the corresponding blockade of each feature. Features 
within the range of level 1 blockade of all amino acids have higher importance value. c, The receiver operating 
characteristic curve (ROC) of random forest model for training, testing and independent validation dataset of L1 
signals of all 20 amino acides. The area under curve (AUC) of different dataset are shown with coloured text label. 
d, Confusion matrix of amino acid identification using RF model. For amino acid Cys (C), and Lys (K), the 
percentage of signals was much lower than 100%, indicating they were identified incorrectly as other amino acids. 
For amino acid Asp (D), and Phe (F), the percentage of signals was much higher than 100%, indicating some other 


amino acids were incorrectly identified as these amino acids. 
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Figure 4. Peptide discrimination by identifying amino acid hydrolysates. a, Schematic illustration of peptide 
hydrolysis using carboxypeptidase Al. Carboxypeptidase Al releases amino acids (except R, K and P) from 
C-terminus of peptide. Then the released amino acids are detected and identified. b, Current trace of real-time 
detection of hydrolysates from peptide EAFNL. ec, The distribution of relative abundance of the blockade of amino 
acid signals identified by machine learning algorithm from the hydrolysate of peptides. d, The identified 
compositions of three peptides. Number in parentheses represents the count of recognized signals. 


